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1.  INTRODUCTION.  The  estimation  of  a  regression  function 
which  is  defined  on  an  interval  of  the  real  line  and  is 
nondecreasing  (nonincreasing)  but  is  not  assumed  to  be  of 
a  particular  functional  form  has  been  considered  in  the 
literature.  Typically,  an  estimate  is  chosen  which  minimizes 
a  particular  objective  function  subject  to  the  appropriate 
montonicity  constraints.  Brunk  (1958,  1970)  considered  the 
weighted  least  squares  estimate  and,  in  the  later 

reference,  demonstrated  its  consistency  and  obtained  its 
asymptotic  distribution.  The  large  sample  distribution  results 
have  been  extended  by  Leurgans  (1979)  and  Wright  (1981).  Rates 
of  convergence  for  this  estimator  have  been  studied  by  Makowski 
(1973)  and  Hanson,  Pledger  and  Wright  (1973). 

The  least  absolute  deviations  (^1)  estimate  was  intro¬ 
duced  by  Robertson  and  Waltman  (1968).  Cryer,  et  al.  (1972) 
also  considered  this  estimator,  showing  its  consistency,  studying 
its  rate  of  convergence  and  comparing  it  to  the  estimator 
by  Monte  Carlo  techniques.  Casady  and  Cryer  (1976)  have  shown 

that  based  on  r  observations ,  the  almost  sure  rate  of  convergence 

-1/4 

for  the  estimator  is  of  order  no  larger  than  r  ,  assuming 

the  underlying  regression  function  satisfies  a  first  order 

Lipschitz  condition.  (This  is  the  same  rate  obtained  by 

Makowski  (1973)  for  the  estimator.)  In  this  paper,  the 

asymptotic  distribution  of  the  l .  estimator  is  obtained  and 

assuming  the  regression  function  has  a  positive  slope  at  a  point, 

the  rate  of  convergence  at  that  point  is  seen  to  be  of  order 
-1/3 

r  .  The  techniques  presented  here  also  apply  to  the  weighted 
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l ^  estimate  provided  the  weighting  function  is  positive, 
continuous,  bounded  and  bounded  away  from  zero.  As  was 
seen  in  the  least  squares  case  (cf.  Wright  (1981)),  such 
weights  do  not  affect  the  limiting  distribution  and  so  we 
only  consider  the  case  of  equal  weights.  (For  a  description 
of  the  weighted  estimate  see  Robertson  and  Wright  (1975)-) 

Leurgans  (1979)  has  obtained  the  asymptotic  distribution 
of  another  estimator  which  might  be  appropriate  if  the  errors 
have  heavier  tails  than  the  normal  distribution.  The  estimate 
she  has  considered  is  defined  to  be  the  slope  of  the  greatest 
convex  minorant  of  a  process  determined  by  smoothly  weighted 
linear  combinations  of  order  statistics.  For  Leurgans* 
estimator,  the  estimator  and  the  l ^  estimator  the  order  of 
the  rate  of  convergence  is  the  same  and  so  the  large  sample 
relative  efficiencies  are  determined  by  the  multiplicative 
constants.  These  comparisons  are  not  discussed  here  since  they 
are  the  same  as  those  for  the  ordinary  one  sample  location 
problem. 


2.  ASYMPTOTIC  DISTRIBUTION  OP  THE  ESTIMATOR.  Let  p  e  (0,1) 

and  for  each  x  e  I,  an  interval  of  real  numbers,  let  D(x)  be 

a  probability  distribution  with  pth  quantile  0(x).  For  each 

positive  integer  r,  let  xrl  <  xr2  <  . . .  <  xrr  be  points  in  I 

and  let  Y  , ,Y  oJ...,Y  „  be  independent  random  variables  with 
rl  r2  rr 

Y^^  distributed  as  •  The  xrl{  are  observation  points  and 

the  Y  ,  are  observations.  (For  the  results  given  here  the 
rk 

number  of  distinct  observations  must  grow  at  least  like  some 
positive  constant  times  r.  If  the  number  of  distinct  observation 
points  is  bounded,  the  work  of  Robertson  and  Waltman  (1968) 
gives  the  asymptotic  distribution.)  The  estimator  proposed 
by  Robertson  and  Waltman  is  defined  at  the  observation  points 
and  any  monotone  extension  to  I  might  be  appropriate.  One 
such  estimator  is  given  by 

(1)  0rU)  =  max  "In  ,  V[xrs-xrt]) 

rs—  —  rt 

where  Qr(A)  is  the  pth  sample  quantile,  that  is  the  [rp]th  order 

statistic,  of  the  sample  comprised  of  those  Y  for  which 

x  .  e  A.  This  choice  of  6  is  constant  on  (x  .  . ,x  .)  for 
rk  r  rj-l  rj 

j  =  2,...,r,  however,  examining  the  proofs  we  see  that  the 
large  sample  results  given  here  are  valid  for  any  nondecreasing 

A 

estimator  which  coincides  with  0  at  the  observation  points. 

r 

If  p  *  1/2  the  monotone  l ^  estimator  is  obtained. 

The  result  that  follows  gives  the  asymptotic  distribution 

A 

of  0  (x_)  witli  xrt  in  the  Interior  of  I.  As  in  the  case, 
r  0  0  2 

A 

the  rate  of  convergence  of  the  estimator,  0  (Xq),  depends  on 
the  rate  of  growth  of  the  regression  function  at  x^.  We  assume 
that  for  some  a  and  8,  both  positive, 


/ 


M 
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(2)  |e(x)  -  e(xQ)|  *  e|x  -  x0|a(i  +  °(1^  as  x  -*•  xQ. 

(Some  implications  of  this  assumption  are  discussed  in 
Wright  (1981).)  Since  the  observation  points  may  be  the 
realization  of  a  sequence  of  random  variables  each  of  which 
has  support  in  I,  we  state  the  conditions  on  them  in  terms  of 
their  empirical  distribution  function, 

Fr(x)  =  cardtk:  xrk  <  x}/r. 

We  assume  that  there  is  a  distribution  function  F,  which  is 
continuously  differentiable  in  a  neighborhood  of  Xq  w.' th 
F’ (Xq)  >  0,  for  which 

(3)  supx|Fr(x)  -  F(x) |  =  o(r"1/(2a+1)). 

To  apply  the  usual  techniques  for  sample  quantiles  we  make 
the  following  uniformity  assumption: 

(4)  sup{ |P{Yrk-e(xrk)<x}-p-nx| :  l£k<r,  r*l,2, . . . , |x| <p}=o(p)  as  p-K). 

Of  course,  (4)  is  satisfied  if  the  Yrk  -  9(*rk)  have  a  common 
distribution  function  G,  with  pth  quantile  0  and  G*(0)  =  n. 

THEOREM.  Suppose  that  6  is  nondecreasing  and  satisfies  (2); 
that  the  observation  points  satisfy  (3)»  and  that  the  observations, 
{Yrk}»  are  independent  for  each  r  and  satisfy  (4).  Then 

(5)  {(a+l)(rF*(xo)n2(p(l-p))"1)o‘/e>1/(2a+1)(0r(xo)-e(xo)) 

converges  in  distribution  to  the  slope  at  zero  of  the  greatest 
convex  minorant  of  W(s)  +  |s|  ,  where  W  is  the  two-sided 

Wiener-Levy  process  with  variance  one  per  unit  time. 

Proof.  The  proof  is  similar  to  those  given  by  Prakasa  Rao 

(1969)  and  Brunk  (1970)  except  that  we  must  approximate  the 

sample  quantiles  by  averages  before  their  techniques  can  be  applied. 
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As  in  the  1 2  case,  we  show  that  ©^(Xq)  is  asymptotically 
equivalent  to  an  estimator  based  on  observations  with 
observation  points  in  a  sequence  of  intervals  which  converges 
to  Xq.  We  make  use  of  the  modification  of  the  arguments  due 
to  Prakasa  Rao  (1969)  and  Brunk  (1970)  given  in  Wright  (1981). 
For  an  arbitrary  c  and  r  sufficiently  large  choose 
a^(r) ,  au(r),  8^(r)  and  &u(r)  so  that 

F(x0)-F(xQ-a£(r))  =  F(x0+au(r) )-F(xQ)  =  2cr~1/(2a+1) 

and 

F(Xg)-F(Xg-6£(r))  =  F(xQ+B^(r)  )-F(xq)  =  cr'1'' ('2a+1\ 


Set 


6*  maxxrse (xQ-a^(r )  ,x0]minxrt.e[x0,x0+au(r)  )Qr(  ^rs^rt^  ' 

/\ 

The  first  step  in  the  proof  is  to  show  that  G^Cxg)  and  0 


are  asymptotically  equivalent  by  showing 

(6)  llmc^«»  lim  suPr+coP{ Vx0)?*9r}  =  °* 

Since  the  pth  sample  quantile  is  a  Cauchy  mean  value  function 

(cf.  Robertson  and  Wright  (1975)),  it  satisfies  the  averaging 

property  used  in  the  proof  of  the  lemma  in  Wright  (1981). 

That  proof  can  be  modified  to  obtain  (6),  but  the  following 

should  be  noted:  with  v  (A)  =  card{k:  x  .  e  A},  I.  the 

r  rk  A 

f  y»  \ 

indicator  of  the  event  A  and  2^  the  sum  over  those  k  for 
which  x  k  e  A, 


{mlny>X0  Qr((X0“B*(r)»y3) 


<  e(x0-e£(r))} c 


U  (Av^( (xQ-6^(r) ,y])>p-(vr( (x0-e^(r)  ,xQ]) 


-1 


-  r),y]W<V8£(r,>^». 


^  * 
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where  AvJ,(iO  =  I<r)  (I<Yrk<e(x0-6^(r>  )  rvrk)/vrt k)  and 


prk  =  p^Yrk-e^xO-0£^r^^*  Furthermore,  if  xrk  >  xQ, 

yrk  -  p^Yrk~0 ^xrk^— 6 *  which  using  (4), 

can  be  bounded  by  p  +  n(0  (x^-fi^r) )-0(x q  ))(1+o(1))  where 
o(l)  represents  quantity  that  does  not  depend  on  k  and 
converges  to  zero  as  r  +  »,  If  xrk  <_  xQ,  then 

urklP+n(0(xo-6^/r) )-e(xrk) )+(e(x0-e^(r) )-e(xQ) )o(l) 

and  so  I((x^B£(r),y]»irk/vr((xO"B£(r)»y))-p+(e(xO-8£(r))“e(xO))o(1)  + 


'(xJ-8£(r),x0](e(x0-^(r))“e(xrk))/vr((x0-B£(r)>x0])- 


Next  we  approximate  the  pth  quantiles  by  averages.  Then, 
since  the  max-min  operation  on  averages  is  the  same  as  the 
slogcom  (slope  of  the  greatest  convex  minorant)  of  a  cumulative  sum 


process,  we  only  need  to  make  slight  modifications  on  the  arguments 
in  the  literature  for  such  slogcom' s.  For  arbitrary  e  with 
0  <  e  <  c,  let  y^(r)  and  yu(r)  be  defined  by  Yu(r)  = 

min{xrk :  F(xrk)-F(xQ)  _>  er-1/^2a+1^  }-xQ  and  Y^Cr)  ■  xQ 

-  max{xrk:  F(xQ)-F(xrk)  _>  er”1^^20141^ } .  Note  that 


(7)  maxxrg e ( xQ-a^ ( r )  .Xg-Y^r)  ]mlnxrte [xQ ,xQ+au(r)  )Qr(  ^xrs  ,xrt^-6] 


<max 


xrse(xO-a^(r)»xO]minxrteCVYu(r)*xO+au(r))Qr(i:XrS,Xrt:i) 


The  large  sample  distribution  for  the  lower  and  upper  bounds  in 
(7)  depend  on  e  and  as  e  +  0  they  approach  the  same  limit. 


The  discussion  for  the  upper  and  lower  bounds  are  similar  and 
so  we  only  give  the  former.  Denote  the  upper  bound  in  (7)  by 
0^.  To  approximate  the  quantiles  by  averages,  we  show  that 


ra/(2a+l)max 


xrse(x0-a£(r)jX0]»XrteCx0+Yu(r)jX0+au(r)) 


n(Q([xrs,xrt])-6(x0))  -  AvJ([xrs,xrt3)|  -  0, 


where  Av^(A)  =  Xj(P-I{y  <e(x  )})/vr(A)-  We  aPPeal  to  the 

r  k —  0 

following  lemma  which  is  a  generalization  of  a  result  given 
in  Ghosh  (1971) : 

Lemma.  Consider  arrays  of  random  variables  (A^:  j=l>2,..., 
Jn,k*1,2,. . . ,Kn»n»l,2». • and  {BJR:  j=l,2, . . . ,Jn,k=l,2, . . . , 
Kn,n=l,2,. . . } .  If  maXjk(Ajk|  is  tight  and  for  each  6  >  0  and 
real  x 


(8) 

P{a"  < x  and  B1?.  >x+6 
jk-  Jk- 

for 

some 

j 

and 

k> 

♦  0 

and 

(9) 

P(Ajk>x+« 

and  Bj k£x 

for 

some 

j 

and 

k) 

-+  0 

as  n 

<»,  then 

maxJk|Ajk- 

•B"k' 

S  0 

• 

Proof.  Because  of  the  tightness  assumption  it  suffices  to 
consider  for  arbitrary,  positive  t  and  M 

P (  I  A j ic— B j ic  I  >  2t  and  _M  <  Ajk  <  M  for  SOme  ^  and 

_<  ^P  { -M+  ( £-1 )  t  < Aj  k£-M+£x , Bj  k^-M+  ( £+1 )  t  for  some  j  and  k) 

+Tp_  P{-M+(£-1)t</C1.  <-M+£t,b”  <-M+(^-2)t  for  some  j  and  k) 
where  L  =  [2M/t]  +  1.  The  above  sums  can  be  made  as  small  as 


desired  for  fixed  M  and  t. 
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Set  nlp  =  vr((-oo,x0-a£(r)]),n2r=vr((-»,x0]),n3r* 
vr((-»,x0+Yu(r)))  and  n4p  =  vr ( (-<*>, x0+au(r ) ) ) . 


We  apply  the  lemma  with  n  =  r# 


y  _  „  t  —  m  v,  Ar-  a/(2a+l) .  „  /  r 

Kr  4r  “  3r*  Jr  "  2r  “  nlr*  Ajk  Avr.(^ 


r  r»n2r’i+1* 


r>n3r+k 


])  and  B*.  =  nra/(2ot+1)  (Q„(  [x 


Jk 


r'L^r,n2r-j+1,Xr,n3r+k])"e(Xo))' 


The  proofs  that  (8)  and  (9)  hold  In  this  case  are  similar 
and  we  only  give  the  latter.  Fix  x  and  note  that  (Bj^.  x} 


is  contained  in 


a 


[xr,n2r-j+l-xr,n3r+k]I<Yr^e(xO)+x/('lrC‘/(20'+1))) 


>  [(j+k+n3r-n2r)p]} 

and  {B1?.  <  x  and  A^,.  >  x+6}  is  contained  in 


Jk  - 


Jk 


<10)  <IC*P.n2p-J+l,x...n^k]V  i  (xM)(J«+n3p-n2p)p-a/(2o+1>.U 


r  *  3r 


where  Ur£  -  >' ’  *"  ^ 

values  of  the  index  l  in  the  last  sum  0 (xQ-a^(r) )<0 (xr£)£0 (xQ+au(r) ) 

and  applying  (2),  both  0 (xQ)-0 (xQ-a^(r) )  and  0(xQ+au(r) )-0(xQ) 

can  be  written  as  0(2c/F' (xQ)  )0tr“ct//^2ct+1^  (l+o(l) ) .  Hence, 
applying  (4),  EVrt  =  Xr"a/(2ot+1)  +  o(r_0t/(2a+1) )  and  since 

n_  -n„  >er2a^2oi+^  (l+o(l) )  for  r  sufficiently  large  (10)  is 
3r  2r— 

contained  in 

dur,„2r-3«,xr,n3r+«<V-E^>  i  £{r“/(20+1>/21- 


With  r  fixed,  define  ^ 


r,n2r*  r»n 


3r 


](Ur£-EUr£)- 
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Xl,j+1  "  Urn2r-j+l  "  EUrn2r-j+l  for  i  1  »2 » *  * '  » Jr  ,Xk+l ,  1 
Urn3r+k  -  EUrn3r+k  for  k  =  1>2-*‘>Kr  and  Xij  H  0 

p 

otherwise.  So  for  sufficiently  larger  r,  {B^<x  and 

A^  >x+6  for  some  j  and  k}  is  contained  in 
J  k 

{lnaxi<j<Jr+i,iik<Kr+i^s=i^t=ixs,t-e,Sra/(2a+1)/2}  and  applyins 
Theorem  6  of  Gabriel  (1977)  the  probability  of  the  latter 


event  is  bounded  by 
-2<x/(2a+l), 


( Jr+Kr+n3r-n2r)maxnlr+l<^<n4r  V(Ur£} } * 


However,  Jr+Kr+n3r~n2r  =  ^cr2ot/^2a+1^  (l+o(l) ) ,  VdJ^,)  £  lEU^, |, 


JvV  -  |l,ur l 


and  for  n^r  +  1  _<  £  _<  n^r,  EU^  was  shown  to  be  of  the  form 
xr-a/(2a+l)  +  o^r-a/(2a+l) ^  Hence,  (8)  holds.  (The  proof 
of  (9)  is  similar.)  Now  we  must  show  that  max.  ,  |  Al\  |  is 

J  k  j  K 

tight.  Since  | El r v  fi,  x  •,  |  is  uniformly  bounded  for 

ilrl-KX  0n 

n^r  +  1  £  t  n^r  and  r  =  1,2,...,  it  suffices  to  show  that 


j  ,k'  Hx 


rn2r-j+l>xrn3r+^<I^rti9(x0>)'EI(Yriie(j'0)))/<'i+k+n3^r') 


I 


is  tight.  But  this  follows  from  an  argument  like  the  one  above 
which  uses  the  result  due  to  Gabriel  (1977). 

Now  we  consider  the  large  sample  distribution  of  0^  = 

n  maXxrse(x0-ae(r)  ,x0-*minxrte|-X0+Yu(r)  ,x0+au(r) ) AVr(  txrs  ,Xrt  ^  ) 
Since  the  proof  is  similar  to  that  given  in  Wright  (1981)  we 

use  the  same  notation.  Let  y  ,  <  y^2  <  ...  <  y  ^  be  the  distinct 

observation  points  in  (xQ-a^(r),  xQ+au(r));  n(r,k)  the  number 
of  observations  at  yrR;  Y=l£=1n( r ,1 ) ;  XrR  =  Av^({yrR}); 
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tr0  =  0,  trk  =  2cD^£=1n(r,£)/y  for  k=l,2,...,A;  D=2(a(xQ)B)~2; 

a2(xQ)  =  p(l-p)/n2;  and  B  =  { (a+1)  (P’ (xQ)  )a/6a2a(x0)  }1/(2ot+1) . 
Define  a  process  on  [0,2cD]  by  Ur(0)  =  0, 

Ur(trk}  =  2cDl£=1n(r,f)Xr£/Y  for  k  =  1,2,..., A, 


and  linear  interpolation  between  the  points  t  Let  j(r)  and 

Ur)  satisfy  yrJ(r)  <  x0  <  yrJ(r)+1  and  yrt(r)  =  x„+Yu<r). 

In  Chapter  1  of  Barlow,  et  al.  (1972)  the  relationship  of 
the  max-min  operator  to  the  left  hand  slope  of  the  greatest 
convex  minorant  of  the  cumulative  sum  process  is  discussed. 


Using  this  relationship  and  denoting  the  slope  from  the  left 
at  x  of  the  greatest  convex  minorant  of  the  graph  of  X(s)  for 
s  e  S  by  slogcom  (x) { (s,X(s) ) :  s  e  S),  we  see  that 
(11)  n-1slogcom(trj ^ ) { (t ,Ur(t ) ) :  0<t<2cD) 

1  1  n-1slogcom(tr£^r)){(t,Ur(t)) :  0<t<2cD) . 

The  argument  starting  on  p.  of  Wright  ( 1 9 8 1 )  shows  that 

ra/(2a+l)Bri""lslogcom(tr,  (r)  ){(tjUr(t))  :  0  <_  t  £  2cD) 
converges  weakly  to 

V  =  slogcom(O) { (s ,W(s )  +  |s|a+1):  -cD  £  s  £  cD}. 


However,  it  should  be  noted  that  since  sup^<£<^  [EI^  <9(x  )} 

-  P  -  n(6(x0)  -  e(xr^))|  *  o(r~a/(2a+1)) ,  rot/(2ot+1)Bn“1EUr(t) 


can  be  written  as  the  sum  of  fr(i)  and  a  function  which 
converges  to  zero  uniformly  in  t  c  [0,2cD]  as  r  -*■  °°. 

In  considering  the  upper  bound  in  (11),  we  first  note 
that  maxi<£<x (yr£~yr£-l^  *  o(r-1//^2ct+1^ )  which  implies  that 
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vr((x0,x0+yu(r)])  =  er2a/(2ct+1)(l+o(l)).  Since  in  the 
arguments  above  we  have  seen  that  trj(r)  +  cD>  ^r^Cr)  cD 
+  eD/2.  We  appeal  to  the  same  weak  convergence  result  as  before 


to  show  that 


a/(2o+l)n  -1 


Bn“  slogcom(tr^^r^ )  { (t  ,Ur(t) )  :  0  <_  t  £  2cD} 


converges  weakly  to 


V(e)  =  slogcom( eD/2 ) { ( s ,W( s )+ | s | a  ):  -cD  £  s  £  cD)  , 


but  we  must  show  that  with  probability  one  the  convex  minorant 

of  W(s)  +  |s|a+1  on  [eD,  cD]  has  unique  slope  at  eD/2. 

However,  if  the  convex  minorant  has  different  left  and  right 

0^ 

slopes  at  a  fixed  point  t,  then  it  agrees  with  W(s)  +  |s| 

at  the  point  t.  This  implies  that  for  h  >  0,  (W(t+h)  +  |t+h|a  +  '1 

-  W(t)  -  |t|a)/h  is  bounded  below  by  the  right  hand  slope  of 

the  convex  minorant  at  t,  but  with  probability  one, 

lim  inf.  „(W(t+h)  -  W(t))/h  =  -°°.  So  with  probability  one 
rr»-U 

the  convex  minorant  has  unique  slope  at  t.  Again,  using  the 

same  techniques  V(t)  converges  weakly  to  V  as  e  -*  0. 

We  now  show  that  V  =  r°‘//<'2a+1^B(6*-6  (xQ) )  converges 

weakly  to  V.  Let  r^  be  a  subsequence  and  V'  a  generalized 

random  variable,  that  is  one  that  may  be  infinite  with  positive 

probability,  with  P{V  <x}-»-  {V’<x}  for  each  real  x  at  which 

rk 

P{V'<x}  is  continuous.  Since  A^Ce)  =  ra//(2a+l)ggH  is  bounded 

above  and  below  by  sequences  which  converge  weakly,  A  (e)  is 

tight  for  each  e  >  0.  So  for  each  c  >  0,  there  is  a  further 

subsequence  r.  , , ,,  with  A  (e)  converging  weakly,  say,  to 

kU)  rk(j) 

A(e).  Let  x  be  a  real  number  in  the  continuity  sets  of  V,V',V(e) 
and  f\(z)  for  a  sequence  of  e  converging  to  zero.  Since  A^(e) 
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and  B  (e)  =  ra// ^2ot+1^B( 0^,-0 (xQ ) )  are  asymptotically  equivalent 
for  each  e  >  0,  for  real  x  and  e  in  the  sequence  mentioned 


above  we  apply  (7)  and  (11)  to  obtain 


P{V' <x) >lim .  P{B 
-  -  j->»  r 


k(  j  ) 


_<x}  =  P{A(e)<x}  _>  P{V(e)<x}. 


Letting  e  -*■  0,  we  see  that  P-fV'^x}  _>  P{V<x}  for  x  in  a  dense 
set.  The  reverse  inequality  can  be  obtained  by  considering 
similar  arguments  for  the  lower  bound  on  V  .  So  V  =  V'  in 
distribution  and  this  part  of  the  proof  is  completed. 

The  only  step  remaining  is  to  show  that  slogcom 
(0 ) { ( s ,W(s )+ | s |  ):  -cD  <_  s  <_  cD}  converges  to  slogcom 

(0)  { (s,W(s)+  j  s  | a+1)  ;  -oo  <  s  <  «},  but  this  is  established  in 
Wright  (1981). 


Casady  and  Cryer  (1976)  have  shown  that  with  probability 
one,  lim  sup^Jr/loglog  r)1/lj(0r(xo)-0(xo))  <  k  <  -  provided  0 
is  Lipschitz  of  order  1.  It  would  be  interesting  to  determine 
if  the  exponent  could  be  increased  to  1/3  when  0’(Xg)  >  0. 

Of  course,  in  this  case,  it  can  not  be  larger  than  1/3. 
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